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Abstract 



We caicuiate the single-particie Green's function for the tight-binding band 
structure, £p = —2tcosp x — 2tcosp y — fi, with a function of chemical poten- 
tial \x for square-lattice system. The form of the single-particle self-energy, 
E), is determined by the density-density correlation function, x(Qi u )} 
which develops two peaks for ji > — 2.5i unlike parabolic band case. Near 
half filling xiQ,^) becomes independent of u, one dimensional behavior, at 
intermediate values of u which leads to one dimensional behavior in T,(p, E). 
However [i < — O.lt there is no influence on the Fermi Liquid dependences 
from SDW instability. The strong p and E dependence of the off-shell self- 
energy, T,(p,E), found earlier for the parabolic band is recovered for fi < —t 
but deviations from this develop for /i > —O.lt. The resonance peak width 
of the spectral function, A(p, E) has linear dependence in ^ due to the E 
dependence of the imaginary part of T,(p,E). We point out that an accurate 
detailed form for E) would be very difficult to recover from ARPES data 
for the spectral density. 
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I. INTRODUCTION 



The Fermi liquid(FL) theory has been used to describe metallic phase in which quasi- 
particle concept develops from a free electron model to study the metallic behavior [|TJ. 
Since the discovery of high temperature superconductivity in the cuprates ||, the normal 
state properties of the quasi-two dimensional cuprates have been investigated experimentally 
and theoretically. The experimental observations such as optical conductivity ||, electrical 
resistivity J4|-|| and angle resolved photoemission spectroscopy(ARPES) |7||| do not follow 
the conventional FL behaviors. The absence of those FL dependences has lead to a proposal 
that FL model for these quasi-low dimensional conductors should be discarded in favor of a 
Luttinger liquid-like model || in which there is a separation of charge and spin degrees of 
freedom in the elementary excitation. 

However Castellani, Di Castro and Metzner have shown that the FL regime is recovered 
for dimensions greater than one, having investigated instability of Luttinger liquid(LL) by 
analyzing correlated fermions with anisotropic hopping amplitudes in one dimensional sys- 
tem |lOj and dimensional crossover FL to LL using analytic continuation for non-integer 



dimensions ||11||. Independent calculations show that arbitrarily small transverse hopping 



kills off LL |l2|,|T3j. They also found no breakdown of the perturbation theory. 

Given that the expected FL dependences are based on parabolic band (a spherical Fermi 
surface model) and a quasi-particle approximation and assuming that the ground state is 
analytically continuous to the non interacting ground state as a function of interaction 
strength, the discrepancy between the FL picture and the experiment has two possibilities. 
Firstly the conventional FL dependences are present but they are limited to anomalously 
low energy and temperature scales induced by interactions among the fermions or by the 
band structure of the non-interacting fermions. In this picture the spectral density of the 
single-particle Green's function is still characterized by a single-particle like sharp resonance 
for momenta very close to the Fermi surface whose width becomes progressively broader 
as the quasi-particle momenta move further from the Fermi surface where the quasiparticle 
approximation breaks down. The second possibility can be viewed as an extreme case of 
the first possibility in which interactions are so strong that the quasi-particle approximation 
breaks down over at least some region of the Fermi surface. In this case the ground state 
is still a FL in that it has developed analytically out of the non-interacting ground state. 



This picture is supported by ARPES measurements [JT4,15| on the underdoped cuprates 



although the quasi-particle peaks remain anomalously broad ever at the chemical potential 
|16|j . In these experiments the Fermi surface seems to develop around a particular direction, 
signaled by a quasi-particle peak in the measured spectral density at low doping. The region 
in which there is a Fermi surface defined by the vanishing of a quasi-particle peak grows 
with doping until the whole surface is established. || The absence of a quasiparticle peak is 
referred to as a pseudogap. Previously this type of data for underdoped cuprates had been 
described in terms of hole pockets which were shown to arise from strong correlation effects 
in calculations of the t-J model |P7^pU||. However a search for these pockets in the recent 
ARPES results || was unsuccessful. 

The leading 2D FL behavior for on-shell self-energy with an isotropic band structure at 

zero temperature, ^ = - is 
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£'(p, fr) = -a& - - 7ie P H 5/2 + (la) 
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£V, = /3£| in lep/eoi + (ib) 

where the cut-off energy £o is isotropic | 21[| . These leading on-shell dependences of imaginary 



part of self-energy had been determined previously by a number of authors . For low- 

density system £ goes to zero so that the region of the generic Fermi liquid behavior is 
limited. The energy dependence in the real part is a mirror of the imaginary part of 

the self-energy, In through Kramers-Kronig relation which gives a T 2 contribution to 
the specific heat in addition to that from zero sound fl2"lf . The |^| 5//2 term is the leading 
zero-sound contribution. This term has a corresponding (£,p — £,th) 3 ^ 2 term in where 
ith = (vzs ~ v f)Pf> determined by the velocity of the zero-sound mode. 2D calculation shows 
no evidence of a breakdown in FL theory and is similar to 3D FL considering a parabolic 
band. The cut-off energy, £ , has a corresponding cut-off temperature, T , which restricts 
the leading correction to FL behavior in E"(p, £ P ,T) so that T behavior will take over 
the FL behavior at low density limit [31]. These cut-offs are determined by quasiparticle 
interactions and the characteristic energy of the non- interacting system, Ej. As the system 
goes to quasi ID band which can be modeled by considering different ratio of t x and t y , for 
example quasi ID organic materials such as (TMTSE) 2 X or BEDT — TTF [P5], van Hove 



singularity moves to the bottom of the band so that particle-hole (ph) pair contribution is 
more important than particle-particle (pp) contribution for low density limit. In this case 
the nature of the qlD FL behavior remains similar to that of 2D FL |26 |. For parabolic 



band structure pp channel which leads Cooper instability is more important than ph channel 
in 3D and vice versa in ID due to the density of states. No particular set of diagrams is 
significant in 2D |27[| , but at least the second order calculation is shared by both of them so 



that functional dependence in either one will represent the FL behavior. The pp channel is 
discussed later in this paper. 

The cuprates have an anisotropic band-structure which can lead to strong band structure 
features in the density of states. This anisotropy of the band structure is used in a number of 
different models to explain the lack of FL dependences in experiment PE|,Z5f|. The anisotropy 
has been invoked to describe transport data in the hot spot |30| and cold spot models. 



In the hot spot model it is proximity to the SDW instability that is important whereas in the 
cold spot model it is proximity to the superconducting phase through pairing fluctuations 
which leads to the anisotropy. More recently renormalization group calculation have shown 
that if the Fermi surface is dominated by a set of saddle points in the density of states 
Landau FL theory, the quasi-particle approximation to FL theory, can breakdown |32| , |33 l. 



The Hubbard model has been employed extensively in theoretical investigations of the 
normal state properties of the cuprates because the tight binding band structure and short- 
range Coulomb correlations result in an effective magnetic Hamiltonian for low doping which 
is consistent with the data. We investigate the FL characteristics of this model as a function 
of chemical potential, fi, from fi = —3t where the Fermi surface is a circle to fi — — O.lt 
, near half-filling, where it is almost a square. A 2D nearest neighbor tight-binding band 
structure, £^ = —It cosp x —2t cosp y — /i, will be used to determine the momentum and energy 
dependence of single properties in FL behavior where t is the nearest neighbor hopping 
constant. Near the half-filling curved M points are connected by flat Fermi surface which is 
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a good model for the band structure of BSCCO in which there is a strongly nested vector 
|34jl . The Hamiltonian includes a repulsive contact interaction which we treat in the weak 



coupling approximation. By introducing a cut-off, q c , for the interaction in momentum 
space the influence of long wavelength interactions on E(p, E) in the ph channel has been 
investigated in this model and shown to determine the functional form of E(p, E) at low 
energies as in the parabolic band case. The evolution of the leading behaviors in real and 
imaginary part is determined by the changes in the bare ph propagator as the Fermi surface 
becomes more anisotropic. First we calculate the contribution to E(p, the on-shell self- 
energy, to the second order in V which gives the functional form of the corrections and 
then consider the effect of the repeated scattering from ph channel which would indicate the 
influence of any instabilities in these behaviors. 

We next calculate the off-shell self-energy, H(p,E), which is probed directly in angle 
resolved photoemission through the spectral density A(p, E). It is found that energy depen- 
dence of the imaginary part of Y>(p,E) leads to a width of the quasi-particle resonance in 
A(p,E) which is linear in £p even when the on-shell T,"(p,E) has the expected £^ln£^ de- 
pendence. We point to the difficulty of extracting the functional form of off-shell self-energy 
from ARPES data on the spectral function. 



II. CALCULATIONS 

From the Hamiltonian 

H = X^J> c p>+ £ Vi^faCp^-frfCiW* ( 2 ) 

p,a p,q,cr,crl 

where £p is the electronic band structure and V is the electron-electron interaction which is 
treated as a perturbation term, the single-particle self-energy, E(p, E) = E'(p, E)+iE"(p, E), 
is given by 

S(p, iE n ) = —T ^2 G {P- Q, lE n - iuJi)V ef f(q, u{) (3) 

where G(p, iE n ) is the unperturbed temperature Green's function and uji is Bose Matsubara 
frequencies. The V e ff includes all the ph diagrams in principle. As a function of filling 
the electronic band structure develops an anisotropic Fermi surface which will end up with 
quasi-one dimensional square Fermi surface at \i — 0. As the system goes to half-filling, 
the density of states at the Fermi surface approaches to the Van Hove singularity, which 
makes ph channel more significant in the low filling system. We are to concentrate on two 
directions in momentum space through out the calculations, (l,0)p and (1, l)p, which are 
extreme directions in the tight-binding band structure. 

The accurate calculation of the density- density correlation function is essential in this 
study. The correlation function, to) = x'(<f, uS) + zx"{q,u), for ph channel is given by 

J U — {^p+q - £p) 

where £p = —2tcosp x — 2tcosp y — \l. The imaginary part of the density-density correlation 
function for the tight-binding band structure is given by 
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via c>o) = ±[dn id, jm2kk£hzM ^ 

AnJ yx J yy | _ 4f / sin 2 f - - sinfe - f ) sin f ] 2 | 



* _ ^-^sin^-fjsinf 
P »~ Sm [ 4tsinf J+ 2 (5bj 

in which subscripts, x and can be switched when q y < q x . The Kramers-Kronig rela- 
tionship is used to calculate the real part of the correlation function. The (1, 0)q or (0, l)q 
direction can be handled by the expression, 

,/ , v 1 Pf(Pl)®(p},max 'Pi) ~ Pf(P2)®(p},max ~ Pi) r^U I U ^2^ f R \ 



where p f (p) = cos x [-f t - cos(p)], p f , max = cos - f t ], p x = sin M sijg^aj ] ~ f and 
p 2 = sin~ [ 4tsi ^^ 2 J + |. Unlike the parabolic band structure the anisotropy of the system 
in the tight-binding band structure restricts the solution of the delta function to different 
areas of the Fermi surface as the chemical potential approaches to half-filling and a direction 
of the system aligns (1, l)q. For a given uo, the solutions of the delta function come from two 
different regions of Fermi surface. Only one of these contribute for the q = (1, 0)g direction, 
whereas both of them survive in q = (1, l)q direction in the long wavelength limit. 

On the ph channel the short wavelengths do not affect the logarithmic behavior in 
Du ^ do contribute to the cut-off energy, £ . At half-filling with q = 2pf, the 
correlation function diverges in both real and imaginary part which will reflect SDW insta- 
bility P5] , |5B| . Hence by restricting a range of the interaction we tune the cut-off momentum, 



q c , in the self-energy calculation to identify the values of q which are responsible for FL 
energy dependence. 

First we discuss results of the 2nd order calculation in the next section which will show 
the FL energy dependences in the function of band structure and we evaluate and discuss 
the repeated scattering case and Stoner instability. Except in ji — > limit, there is no 
nesting phenomena from short wave length which solely contribute the cut off energy, £o 
without affecting the log behavior. We consider two different regions of fi separately as in 
anisotropic 2D FL and quasi ID FL. 



A. 2nd order calculation of £p) 

The self-energy at zero temperature for second order calculation is given by 

£(£ E > 0) = £ e(E - ^)e{^)V q 2 X {q, E - ^) (7) 

where x(<L w ) is the density-density correlation function or response function and V$ = 
VQ(q c — \ q\) in which q c = tc and V = t throughout the calculations. The £| In behavior in 
the imaginary part of the self-energy comes from the linear energy dependence in x"(q,u)). 
In the self-energy calculation the energy dependence is mainly determined by the frequency 
dependence of the density-density correlation function rather than the step functions in 
the expression for the self-energy. Since the second order contribution is common for both 
of ph and pp channel, the second order calculation reveals the functional dependence of FL 
behavior as a function of chemical potential which we break in two regions, 2D and quasi-lD. 
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1. Anisotropic 2D Fermi Liquid : —At < jj, < — t 

As the chemical potential, /j,, goes to zero, flatness of Fermi surface along (1, l)p direction 
develops. Although the Fermi surface develops flat, quasi-one dimensional, regions for fi ~ 
— t, no signature of quasi-one dimensional behavior is seen in S(p, E). The figure |2| shows the 
evolution of the response function with a function of \i in the directions of (1, 0)q, (1, 2)q and 
(1, l)q with q = 10~ 2 7r. The cut-off energies, u* and associated with the contributions 
given momentum q are from the two different Fermi surface contributions through the delta 
function for a given momentum q. Those contributions give two different locations of the 
divergent peaks, u* and ct>|. As q points in the (1, 1) direction, uj\ approaches zero frequency 
so that the FL linear behavior is limited to small energies. However even in (1,1)? direction, 
the linear piece survives and H(p,E) has FL behavior for \i = —t. As uj{ goes to zero 
energy, the linear behavior in x"(9> a; ) is limited to a small u region above which x"( < fi a; ) 
is independent of u as fi — > as in one dimensional system. This will be discussed later. 
Figure |2] shows the anisotropic band structure effects, oj* and ui^ feature, in particle-hole 
correlation function which first appear around \i = 2.5t. Except fi ~ case, the x"^^) 
follows u>/\q\ at small u which leads to ^ In ^ form for T,"(p,^p). 

The log behavior of the imaginary part of the on-shell self-energy comes from long wave 
length calculation. The figure ^ demonstrates contributions \q\ < in x(g, cj) to the 
imaginary part of on-shell self-energy with six different sets of \q\ values given by 10 _ra 7r < 
\q\ < 10 _( - n_1 V where n = 6,5,4,3,2,1. For small values of \^/t\ the contributions from 
longer wavelengths provide the ^ln^ dependences whereas shorter wavelengths lead to 
contributions to the £| dependence. The log behavior is restricted by small energy and the 
\q\s which are long wavelength. The anisotropic 2D systems with —At < // < — t shows 
generic FL behavior with increasing anisotropy in p. The correction term is enhanced by 
tuning the chemical potential close to zero and will be discussed in the next section. The 
cut-off, £o; and f3 for the two directions of p reflect the growing anisotropy of the Fermi 
surface as changes from — 3t to —t. 



2. Quasi- ID behavior : —t<fi<0 

The purely one dimensional response function with a tight-binding band structure, £ p = 
— 2t cosp — fi, is 

! efr, - |f -sin-^^D-e^- if + sin- 1 ( i ^ I )|) u 

X ^ U) = 2 |4tsm V l- W | 6(1 - W ) 

(8) 

where pf = cos~ 1 (—fi/2t) is the Fermi level for one dimensional system. In the small 
frequency limit there is no energy dependence in x'Id which has corresponding logarithmic 
divergence in x'idi 

^-^ZS^jM^gl. (9) 
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At static limit this divergence at q = 2pj, so called perfect nesting vector, drives the Peierls 



instability via Kohn anomaly which leads to the metal-insulator transition [[37 . 

Coming back to 2D problem, at /z = the system is one dimensional in that there is 
perfect nesting between parts on the Fermi surface along (1, l)p direction. Nesting vector, 
q = (1, l)ir, gives a form 

x "{ir, n,Lu ;t i = 0) = ^sign{u)N{u/2) (10) 

where N(u) is density of state which has logarithmic divergence at u = 0, which breaks 
quasi-particle picture. For ji = the constant energy dependence at long wavelengths in 
x"(q,uj) which is one dimensional effect lasts up to zero energy. Using the Kramers-Kronig 
relationship this energy dependence for small \q\ in x"{Qi u ) leads to a lncu divergences in 
the real part, x'iQi u )i signaling the break down of RPA. The large \q\ behavior is shown in 
figure ^| where it is clear that the In a; divergence in cu) is rapidly destroyed as /i goes 
away from half-filling. 

For \q\ < 10 _1 7r structure of x" ■> two divergent peaks identified by uj\ and u^, follows 
the discussion in previous section. Due to rapid development of frequency independence 
in x" f° r l° w energies as fi — > 0, the quasi ID corrections to ^ln^ behavior of £p) 

at intermediate values of u have a £^ dependence in the quasi-particle life time. This is 
shown in figure [7] for /i = — O.lt in which the Fermi surface is flat over an extended area. 

3/2 

The correction sets in at £^ ~ 0.02t, which limits generic FL behavior to small region. 
As \i goes away from half filling, the logarithmic divergence in u dies out rapidly as in the 
figure |5| so that the system is far from a Peierls-like instability. The second order calculation 
shows that quasi-particle life time has generic FL dependence with /i ^ with a term 
as a correction to £|ln£p- when \i < 0. 



B. RPA calculation of 

1. Long wavelength contribution 

With /i ^ the long wavelength limit contribution which is responsible for £| In £^ or 

{^J 2 for /i ~ survives unless the strong interaction limit is considered in RPA. The effective 
interaction in equation (2) for repeated scattering involves in two independent channels, the 
symmetric and the antisymmetric channels, followed by spin exchanges of and 1 and is 
given by 

v m ,,i - 1 YMMl i 3 v^x&u) m > 
K//few) -5i-^(,» + 2i-^(,» 

where V s = V q and V a = —V q . In RPA the imaginary part of self-energy, J2"(p,^), is just 
enhanced from the second order calculations as in figure p, — — O.li case includes 
corrections in (l,0)p direction. The real part, ^), has ^|^| correction to £p term. 
In the figure ^| the £p term is extracted and comes from long wavelength limit. The £p|£p| 
term is a mirror of £^ln£^ term in £p) through Kramers-Kronig relationship and the 

deviations from the fits are consistent to both of £^>) and S'(p, ^) graphs. The on-shell 
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calculation of ^1^1 term in (1,0) direction with fi = —O.lt stops at £p ~ 0.05t because the 
phase space runs out for that direction and the correction to the £p|£p| term sets in earlier 
than (1,1) direction. The coefficients of the £p|£p| terms for the two directions are different 
due to anisotropy of the system. 

The figure [It] shows the functional dependence of the cut-off energy in imaginary part 
of the self-energy, £ ,p, and the coefficients, f3p, for the tight-binding band structure. The 
anisotropy from f3p grows as the band is filled, but near half filling the trend changes rapidly 
due to the flatness of the band structure. The cut-off energy shows anisotropy even at 
fj, = —3t which is different from the figure |2] in which the anisotropic effects of splitting the 
two different peaks in x"(<f.u;). 



2. Instabilities 

As discussed above the correlation function itself does not indicate any instability unless 
fi = as in figure [5[ For the repulsive interaction antisymmetric channel in effective inter- 
action could have an instability indicated by a pole at u = in the effective interaction. 
The solution for it is given by 1/V$ = x'(<fj ^ = 0) in RPA in which there is divergence with 
perfect nesting vector at fi = so that there are solutions for all |V^|'s. Within the RPA the 
question is how fast the system goes to ID as a band approaches to half-filling. 

In parabolic band structure, the real part of the correlation function in static limit is 
given by 

""'"^^W^ttS 1 (12a) 



J(q/2p f ) 2 - 1 

X2 D (q) = N 2D (0)[l ~ \ q /2p f ) g «g/ 2 P/) " ^ ( 12b ) 

where N(0)s are the densities of state on Fermi level for each dimensions and 9 is a step 
function. In 2D the correlation function has a cusp and a discontinuous derivative at q = 2pf 
due to the step function whereas 3D case varies smoothly. Since the Fermi surface consists 
of two points in ID, the nesting condition £p — £ P +2 P/ is satisfied over the entire Fermi 
surface resulting in a log divergence at q = 2pf which is responsible for the breakdown of 
the quasiparticle picture leading to a metal-insulator transition. 

The figure [IT] the static correlation function in (1, l)q direction has a developing diver- 
gence at |g| = 2p ( j' 1 ^ as the system changes 2D to ID. The function is repeated at \q\ — 7r due 
to periodic lattice where \q\ = tc is at the middle of the dip right next to the peak. \i = —3t 
which has the Fermi surface close to circle shows constant q dependence with densities of 
state 

N(fi) = lim lim x'(<f, u) (13) 

lj— >U 



s 



as in parabolic band case. The phase diagram, V$ vs. //, is calculated by 1/V$ = 
x'(2pf,2pj,u = 0). With a interaction strength of V q = t, the Stoner instability does 
not occur until /i ~ — 0.03t. In the calculations discussed in this paper V q = t dotted line in 
figure ^ so that the instability of the system is absent except when /i ~ for V q < t. 

This conclusion is supported by a recent renormalization group analysis of the Hubbard 
model by Halboth and Metzner |38| who found that the tendency towards the antiferro- 
magnetic nesting instability is limited to fi ^ — O.lt for U — t as evidenced by a growing 
spin susceptibility. They found in fact that d x 2_ y 2 superconductivity fluctuations provide 
the dominant susceptibility further from half filling but even these were suppressed beyond 
fi ~ — O.Olt. The influence of superconducting fluctuations are missing in the present calcu- 
lation |39| . 



Going beyond this weak-coupling region is difficult because of the uncontrolled nature of 
the RPA. One such attempt is the so-called self-consistent fluctuation exchange approxima- 
tion (FLEX) which however is equally uncontrolled. It has been pointed out by Trembley 
and coworkers that the absence of vertex corrections in the FLEX leads even more severe 
breakdowns of sum rules than does the RPA discussed here. Vilk and Trembley have pro- 
posed an improvement in RPA to ensure that the sum rules are satisfied. In the approach 
the influence of the cross channels in the ph expansion on the vertices is mimicked by letting 
the irreducible vertices in the charge and spin channels be independent of each and chosen 
to satisfy sum rules. In this approach SDW fluctuations above T = lead to the breakdown 
of the quasi-particle approximation for U = 4t. These approaches have been limited to 
finite size systems and so miss the long wavelength interactions responsible for the leading 
corrections of interest in the present work. 



3. Contributions from the Particle- Particle Channel 



Fukuyama et al. have investigated the contribution to E) for the pp channel us- 
ing the parabolic dispersion. The contribution to E) in second order in the interaction 
can be thought of as being in either the ph or pp channel. In the second order calculation 
the difference between the two channels is that £jdn£^ contribution comes from propagating 
particle or hole pairs with momenta q ~ 2pf in pp channel whereas it comes from long- 
wavelength particle-hole pairs in ph channel as we shown above in figure |3|. Taking the 
second order diagram to be in the pp channel contribution, its contribution is screened by 
repeated scattering in that channel. In contrast to the case of an atractive interaction there 
is no Cooper instability but instead the contribution is reduced by the repeated scattering. 
Unlike the situation in the ph channel the result is not simply an enhancement or reduction 
of the contribution from the second order diagram giving £pln£ p contribution to £ p ) 
again. Rather there is a strong momentum dependence in real part of the pp propagator 
at q ~ 2pf which eliminates the £pln£ p dependence. As a result the effect of including the 
pp channel in a calculation of E) is to remove the contribution of the second order 
term in the interaction from leading dependences. The pp contribution to E) in a tight 
binding bandstructure remains to be worked out in detail. However the result is unlikely 
to be quantitatively different from the result for the parabolic band. Indeed for particle- 
pairs/hole-pairs with net momentum Q = (it, tt) it is easy to show that the propagator 
K(Q,ui) behaviors as K(Q,u) ~ ^_p 2 /(//) where f(fi) is a polynomials in fj, and gives a 
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vanishingly small contribution to E(p, E) where the consequences of the almost nested Fermi 
surface would be expected to be greatest. 

There have been a number of investigations of collective modes in the pp channel. En- 
gelbrecht and Randeria j|lj] found a collective mode below the two particle continuum in a 
bound state of hole pairs which arises because of the finite density of states at bottom of the 



E_ 

E a 



5/2 



band in 2D. The contribution of this mode to the self-energy is given by E"(p/, E) ~ fx 
where /i is the chemical potential and E a is an energy scale describling the interaction. So 
this mode does not contribute to the leading dependences in H"(p,E). Yang |^2[ and more 
recently Demler and Zheng H3J have found evidence for a different singlet and triplet col- 
lective mode in the pp channel which exists in a narrow region of momenta near (it, tt). The 
narrow range of momenta suggests that this mode can have little effect on the low energy 
properties of E(p, E). 

In summary neither the continuum nor collective mode contributions from the pp channel 
contribute to the leading dependences in E(p, E). 



C. Off-Shell Calculations, E): Spectral Function 

The spectral function can be measured by ARPES and is given by 

A(P, E) = f , (14) 

KP } [E-^-^ip,E)^ + [E"{p,E)} 2 K } 

where the off-shell self-energies are involved. For large cut-off momentum, q c , and parabolic 
band structure the leading contributions to the imaginary part of off-shell self-energy which 
determines the width of the spectral function has a form |24| 



E) = C[{% + 2&E - &)} \n(max[&, \E\}) + (E - H\E - &)] (15) 

where C is a constant. In the figure [l| [E"(p, E) — H"{^ ! p)]/ (E — £p) vs. In \E — £p\ is plotted 
for E — £p < 0. In the limit (E — £p) — > we find 2C^p-ln |^| consistent with equation (15) 
for given p. For /i = — 2t the Fermi surface is close to parabolic band structure and the 
off-shell self-energy has a (E — £p) independent behavior for small values of (E — £p) as 
in the equation (15). For fi = —O.lt the Ti"(j)f,E) is consistent with the equation (15). 
However for p = 0.95p/ although E"(p, £p) still has the ^ln^ dependence, E"(p, E) is seen 
to deviate from the form given in equation (15). The dependence on E in figure shows 
a £' 2 ln|£'| behavior for fi = —2t only for p = Pf and the more complicated dependence 
shown in equation (15) is found for \p\ ^ pj. This is the case for other values of fi also. 
This indicates that the form E 2 hiE for H(p,E) frequently used in the literature is a poor 
approximation ELM to the more complicated function of p and E. 

ARPES data is analyzed in terms of the single-particle spectral density. The extent to 
which a quasi-particle resonance is well defined is characterized by the width of the resonance 
at half-maximum. This is frequently thought of as a measure of the quasi-particle life-time 
where momentum dependence is taken as a test of the Fermi liquid character of the material 
under investigation. When the quasi-particle approximation is applied, the spectral function 
is 
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= W=WTT* (16) 

where Ep is the location of the resonance and is half-maximum width of the peak. Here 
we apply this analysis to the Hubbard model. 

The spectral functions are calculated for \i = — 2t and fi = —O.lt cases and the widths of 
the resonance at the half of the peak, Tp, are plotted as a function of £p < 0. The spectral 
functions are picked along r to Y point and to M point. The quasi-particle peaks are well 
defined and as p — > pf resonance gets sharper for (1, 0)p and (1, l)p directions in the figure 



15| for n = —It. The fx = —O.lt case in the figure [U] has the same feature except growing 
anisotropic feature for (1, l)p direction in deep inside of the Fermi surface. There is no sign 
for the pseudogap with the parameters used in this calculations for \x = —It and \i = —O.lt. 
As p — > pf the resonance peak does not show symmetry about E = Ep since the imaginary 
part of self energy is always zero for E — 0. A(p, E) is only symmetric about E — for 
\p\ = \pf\. The smallest deviation from pf leads to strongly non-symmetric A(p,E). This 
indicates that one can be misled by assuming symmetry for the ARPES data cutoff by Fermi 
Dirac distribution function to get the single-particle properties near Fermi surface. 

If we expand A(p, E) near Ep — £p+ E) which is approximately the location of 
the resonance peak since £(p, E) is continuous, the width of half maximum is given by the 
on-shell self-energy 

r ggg^Cg E p) / 17 n 

""~ 'l + [ Z0 dE"(p,E)/dE\ E=E ^ 



where zp — [1 — dY!(p, E)/dE\ E=E J l 1 unless E) has a strong energy dependences. In 
this calculation we find 1 — zp < 3% showing that there is very little frequency dependence 



in the real part, Ti'(p,E) ~ Ti'(p,^p). The figure |17] shows the agreement between the half 
maximum width and the on-shell self-energy at low energies, Tp ~ £p), for /i = —2t 
and /i = —O.lt. As in the figure the generic FL logarithmic behavior is restricted to low 
energies. For /i = — 2t case Tp shows linear behavior in £p at higher energies which is not 
extrapolated to zero energy. Clearly the assumption that a deviation from Tp ~ £| In £p or £| 
is an indication of non- Fermi liquid is not correct. We note that no sign of the complicated 
£p and E dependences of T,(p,E), shown in equation (15), comes out of this analysis of 
ARPES data. It seems to us unlikely that it will be possible to recover these dependences 
without prior knowledge of the functional form to which to fit the data. 



III. CONCLUSION AND SUMMARY 

We have demonstrated that anisotropy in band structure does not qualitatively change 
FL dependencies. We have done that by calculating single-particle properties in the Hubbard 
model in which the shape of the Fermi surface is changed from a circle to a square as the 
value of the chemical potential is varied from fi = — 3t to /i = —O.lt. In these calculations 
we have used a weak coupling approximation in which the magnitude of the interaction, V q , 
is far from the values for the SDW instability except for fi > — 0.5t. Numerical coefficients 
in the functional forms of the imaginary part of the self-energy reflect the anisotropy of the 
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band structure but it is not until the Fermi surface is almost a square that the first sign 
of quasi-lD behavior, £^j/ 2 , appears. As p — > the magnitude of £p>) increases with 
with increasing anisotropy. Also, in the p — * limit, an a; independent region in x"{Qi u ) 
grows, which is responsible for t^ 2 behavior in E"(p, £p) in (1, 0)p direction. The ^|^| term 
is also found in Tf(p,£p) as in parabolic band case. The instability of the system is slowly 
developed as the system goes to ID except p ~ 0. 

The off-shell self-energy is shown to have a similar p and E dependence as in the parabolic 
band with deviations as the Fermi surface developes flat regions. There is no evidence for 
the pseudogap phenomenon for the parameters used here, V — t and p < — 0. It. In this 
region we confirmed the half maximum width of the resonance peak for spectral function is 
approxmated by on-shell self-energy. Asymmetry in the resonance peak grows for \p\ — Pf 
due to H"{p,E = 0) = and for deep inside of Fermi surface for p = —O.lt and (1, l)|p| 
direction. In the tight binding band structure quasi-particle picture is well defined for 
p < -O.lt and \V$ < t. 

As one goes further away from this parameter region it becomes difficult to justify the 
simple RPA weak-coupling approximation used here. SDW and superconducting fluctuations 
are likely to play a role as other authors have pointed out. Their effects on the leading 
dependences from long wavelengths in Y,{p,E) remain to be investigated. 
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FIGURES 
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FIG. 1. The band structures for the tight-binding dispersion, £p = —2tcosp x — 2tcosp y — ji 
with different values of ji. As band is filled up, the parabolic Fermi surface is changes toa square 
surface. For fi = — O.lt, there is a large region of flat Fermi surface. 
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FIG. 2. The correlation functions as a function of \i and directions for |<f| = 0.017T. (a) x" as 
a function of fi. Two peaks are appeared in \" f° r the anisotropic band structure, (b) x" as a 
function of directions in q. As q points (1, 0)<f, uj\ gets closed to u)2- (c) The locations of peaks in 
the correlation function, oj\ and u)2- As fj, — > the location of the first peak approaches to zero. 
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FIG. 6. The correlation functions for q = 0.0l7r with different /xs. As \x — > to independent 
behavior in low frequency is develops which is characteristic of ID feature at long wavelengths and 
the system shows ID effect. 
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H=-0.1t 
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FIG. 7. The imaginary part of self-energy in second order calculations with V = t shows £ 3 / 2 
behavior which is coming from frequency independent x" ■ This crossover from £jdn£p- to £jV 2 
appears in (l,0)p direction. The dotted lines are fitted lines. 
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FIG. 8. The imaginary part of self-energies in RPA with V = t for (1, 0)p and (1, l)p directions 
with (a) n = —2t and (b) ji = —O.lt. For \i = —O.lt case, behavior for (1, 0)p direction appears 
as in 2nd order calculation. Both of the directions show the logarithmic behavior. 
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FIG. 9. The real part of self-energy calculated in RPA with V = t shows behavior correc- 
tion to a^p-term. The slope 2 indicates the power of the fit. The deviation from the ip\£,p\ behavior 
comes from the long wavelenght quasi ID frequency dependence of x((f,^)- 
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FIG. 11. The real part of the correlation function in static limit. The system changes 2D to 
ID as \i — > 0. A divergence grows at q = 2pj as \i — > in (1, l)q direction. 
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FIG. 12. The phase diagram of the system from 1 — V q x'(2pf, 2p/, 0) =0. For any strength of 
interaction at fx = the quasi particle picture breaks down due to the divergence. With V q = t, 
there is no SDW instability for fi < — O.lt. 
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FIG. 13. The imaginary part of the off-shell self-energy for E < and p < pf in second order 
calculations. E — > ^ limit T>"(p,E) has constant energy term, £^ln£^, for fi = —2t in ln\E — ^\ 
scale in the case of p ^ pj. There is a significant deviations in E) for n = —O.lt in the case 

of p^p f . 
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FIG. 14. The graph shows that the off-shell self energy follows E 2 InE only when p = Pf even 
for fj, = —2t which is close to parabolic band structure. 
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FIG. 15. The spectral functions in repeated scattering calculations for fi = —2t in (1, l)p (a) 
and (1, 0)p (b) directions. The quasi-particle peakes are well defined as p — > py. Asymmetry in the 
resonance peak apears when p ~ pf. 
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FIG. 16. The spectral functions in repeated scattering calculations for [i = —O.lt in (1, l)p (a) 
and (1, 0)p (b) directions. The quasi-particle peakes are well defined as p — > pj. Asymmetry in the 
resonance peak apears for (1, l)p in deep inside of Fermi surface. 



30 



|U=-2t, (l,l)p 





U.UH-O 




n r\Ar\ 

U.U^fU 




0.035 


■<— > 


U.UoU 






o 




1 ■ 




'3 








u 












s= 
i— ( 


0.015 




U.U1 (J 




0.005 




n nnn 




U.I 4 




0.12 


■(— > 


0.10 


o 




■4— > 




5=1 


0.08 


=3 










0.06 


5= 










0.04 




0.02 




0.00 




-1 



x i I 1 1 s 1 1 1 1 1 

\ 

XX \ 

\X \ 

\^ X 

\ N \ 

V\ \ 
\X \ 


i 1 i 1 i 1 i 1 

— r 

p 

F.L. log-fit 
ImZ(p4 p ) 


— / 

Linear behavior 




- a) 

1,1,1,1,1, 


\\ 
X^> 



H=-0.1t, (l,l)p 



\X 1 l 1 l 1 


i 


1 1 1 1 1 1 1 1 


X 
\ 
X, 
V 

x 
x 

X 




r 

p 


x 
x 




F.L. log-fit 


x 

X 
X\ 




imi:(p,y 


Xx 

x 

\ 




b) 


\ 
\ 

X 




i , i , i , .- -"T" 


1 





-1.0 -0.9 -0.8 -0.7 -0.6 -0.5 -0.4 -0.3 -0.2 -0.1 0.0 



FIG. 17. The comparison between and on-shell self-energy for n = —It (a) and ji = — O.lt 
(b). The FL behavior, £j-ln£p-, is restricted to low energies. 
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